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ABSTRACT 

Data from 58 strong lensing events surveyed by the Sloan Lens ACS Survey is used to esti- 
mate the projected galaxy mass inside their Einstein radii by two independent methods: stellar 
dynamics and strong gravitational lensing. We perform a joint analysis of these two estimates in- 
side a model with three degrees of freedom with respect to the lens density profile, stellar velocity 
anisotropy and a possible line-of-sigh (l.o.s.) mass contamination. We also consider two possibili- 
ties for the lens light distribution (Jaffe and Hernquist profiles). A Bayesian analysis is employed 
to estimate the model parameters, evaluate their significance and compare models. We find that 
the data favor Jaffe's profile over Hernquist's, but that any particular choice between these two 
does not change the qualitative conclusions with respect to the features of the system that we 
investigate. The data do not strongly constrain the l.o.s. mass contamination, yielding no contri- 
bution from secondary matter at maximum likelihood and an average contamination of the order 
of 10%, when considering the full range of probabilities. The density profile is compatible with an 
isothermal, being sightly steeper and having an uncertainty in the logarithmic slope of the order 
of 5%. We identify a degeneracy between the density profile slope and the anisotropy parameter, 
but we encounter no evidence in favor of an anisotropic velocity distribution on average for the 
whole sample. 

Subject headings: dark matter — galaxies: elliptical and lenticular, cD — galaxies: kinematics and 
dynamics — galaxies: structure — galaxies: fundamental parameters — gravitational lensing 



1. Introduction 

The observation of strong gravitational lens- 
ing events has allowed many studies about the 
mass, density profi le and structure of t he galaxies 
that act as lenses (IBolton et al. 2006; Treu et al.l 



2006 



Koopmans et al. 2006c iGavazzi et alj 2007; 
Bolton et alj 120071 l2008aUbl; ICzoske et aTlhoOS ; 
Treu et alT 2008 ). which could also have impor- 



tant implications for dark matter and cosmol- 
ogy studies - see als o Kochanek's contribution in 
Meylan et al. ( 20061 ) for a more general review in 
the field. Therefore it is fundamental to control 
for possible systematic effects, such as line-of-sight 
(l.o.s.) mass contamination. However most of 
these works consider a free path from source to 



lens and from lens to observer, even though there 
are observational and theoretical suggestions that 
the l.o.s. mass contamination may be significant. 
Bar-Kana ( 19961 ) investigated theoretically the 



effect of the Large-Scale Structure on strong lens- 
ing events, f i nding that it can be significant. 



Keeton et al.1 fj 19971 ) observed that external shear 



due to galaxies and clusters associated with the 
primary lens or along the l.o.s. can be an im- 
portant perturbation in individual lens models. 
Premadi fc Martel ( 20041 ) used ray-tracing to in- 



vestigate the effect of density inhomogeneities 
along the l.o.s. of strong lenses and concluded 
that it can be of ord er 10% on the ma gnifica- 



tion of the sources. iDalal et al.l (|2005f ) found 
that significant errors can arise from l.o.s. pro- 
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jections when using giant arcs generated by clus- 
ter strong lensi ng to constrain the cos mologi- 
cal parameters. IWambsganss et aL ( 20051 ) noted 
that secondary matter along the l.o.s. of strong 
lenses can lead in some circumstances to an over- 
estimate of 10-15% of the primary lens mass 
i f ignored. Using spectroscopy, various groups 
Tonrv fc Kochanekl I2OO0I iMomcheva et al.ll2 006: 



Fassnacht et a l. 2006; Willi ams et al.ll2006l : lAuger et al 
20071 ) discovered a sig nificant l.o.s. effect on 



strong lens galaxies. iMoustakas et al.l (J2007) 



observed that even in under-dense local envi- 
ronments, the l.o.s. contamination may give a 
considerable contribution to galaxy-scale strong 
lenses. Usi ng ray-tracing t hough t the Millennium 
Simulation, iHilbert et ail (|2007l) determined that 
strong lensing lines-of-sight are biased towards 
higher than average mean densities, contributing 
a few percent to the tota l surface density, and 
Puchwein fc Hilbertl (|2009h found that secondary 



matter along the l.o.s. has a large effect on the 
strong-lensing optical dep th and the c ross-s ection 
for cluster strong lensing. Treu et al. ( 2008f ) mea- 
sured the over-density of galaxies around SLACS 
lenses and observed that typical contributions 
from external mass distributions are of order of 
few percent , but r eaching 10-20% in some cases. 
Faure et al. ( 20081 ) considered strong and weak 



lensing observations in the COSMOS survey and 
compared with simulations, finding that strong 
lensed images with large angular separatio n were 
i n the densest regions. On the other hand, [Auger 
(|2008l ) did not find an over-density of photomet- 
ric sources along the l.o.s. of a limited sample 
of SLACS strong lenses in comparison with other 
SDSS massive early- type galaxies and interpreted 
that as evidence against a possible l.o.s. contami- 
nation. 

In this paper we use two independent galaxy 
mass estimate methods, strong gravitational lens- 
ing and stellar dynamics, to examine the presence 
of l.o.s. mass contamination in the lens set and 
its effect on the determination of the density pro- 
file. We use all the suitable events in the SLACS 
sample, considering realistic brightness functions 
for the lens-galaxies, and incorporating our prior 
ignorance on their velocity anisotropy. 

In Section [2] we describe how we calculate the 
mass of SLACS lenses using strong gravitational 
lensing and stellar dynamics, and how we deter- 



mine the l.o.s. mass contamination. In Section [3] 
we show our results, which are discussed in Section 
IH together with our conclusions. 

2. Data and Methods 

The analysis in this paper is based on the com- 
parison of galaxy masses calculated through two 
different methods: gravitational lensing and dy- 
namical analysis. In sec. 2.1 we present the data 
used in the analysis, collected from the SLACS 
survey. 

In sections 2.2 and 2.3, we discuss the lens- 
ing and dynamical mass determinations, respec- 
tively. We will assume simple models for the 
galaxy mass distribution (e.g., spherical symme- 
try, power-law density distribution) because they 
have few free parameters and allow to illustrate 
well the two methods. For a similar approach see 
Koopmans et al. ( 2006h . 

We want to examine whether the two mass esti- 
mates are indeed equivalent and/or if there is ev- 
idence of systematic differences between them. In 
particular, we consider the plausibility of a con- 
tamination affecting the lensing mass. Sec. 2.4 
presents a Bayesian analysis of this problem. 

2.1. Data 

The selected set of galaxie s is part of the Sloa n 
Lens ACS Survey, SLACS (|Bolton et al.l 120061) . 
which is a Hubble Space Telescope (HST) Snap- 
shot imaging survey for strong gravitational galac- 
tic lenses. The candidates for the HST imag- 
ing were selected spectroscopically from the SDSS 
database and are a sub-sample of the SDSS Lumi- 
nous Red Galaxy (LRG) sample. 



We u se data comp il ed fro m iKoopma ns et al 




Gavazzi et alJ (|2007l ) and iBolton et al 



, constructing a sample of 58 strong gravi- 
tational lensing events where the lenses are iso- 
lated early-type galaxies (E+S0). Data from 
SLACS is especially suitable for joint strong lens- 
ing and dynamical analysis because they allow 
precise determination (5% error) of the Einstein 
radius for each lens-galaxy in a relatively homo- 
geneous sample of early-type galaxies. And, at 
the same time, SDSS has precise stellar veloc- 
ity dispersion measurements (6% average error) 
for the lenses, as well as redshifts for lenses and 
background sources. 
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For each lens system we are interested in the 
redshift of the background lensed source z s , the 
redshift of the lens zi , the average stellar velocity 
dispersion inside an aperture a ap , the effective an- 
gular radius 9 e f and the Einstein angular radius 
9e- The sample average values for these quanti- 
ties are (zi) = 0.2, (z s ) = 0.6, (<J a p) = 250fcr7is _1 , 
(fi ef ) = 2.2", and (9 E ) = 1.2". 

The source and lens redshifts were determined 
from the SDSS spectra, and the stellar velocity 
dispersion corresponds to the light-weighted aver- 
age inside the 3" diameter SDSS fiber. 

2.2. Lensing Mass 

The estimated projected mass inside the Ein- 
stein radius Re = OeDl, is given by 



Mr = 



c 1 D L D S 2 



4G D LS 



(1) 



where D[ L g L gi is the angular-diameter distance 
of the lens, source, and between lens and source, 
respectively. These distances are calculated as- 
suming a redshift-distance relation derived inside a 
chosen cosmological model that in the present pa- 
per is a concordance ACDM model with H m = 0.3, 
il A = 0.7 and h = 0.7. 

The Einstein radii were determined from HST 
images using strong lensing modeling of the 
lenses and reconstruction of the unlensed s ources 
(|Koopmans et al.ll2006l ; iGavazzi et al.ll2007l ). The 
uncertainties on 9e were reported to be around 
5%, so we use this value for all Einstein radii 
when calculating the error on Me- Note that 
the lensing modeling uses a Singular Isothermal 
Ellipsoid (SIE) mass model, but the resulting pro- 
jected mass distribution is parameterized by an 
Eintein radius so that the enclosed mass in the 
projected ellipse is the same that would be en- 
closed in a projected circle from an equivalent 
Singular Isothermal Sphere. This is the radius 
we adopt here. Indeed, the Einstein radius deter- 
mined this way is a robust attribute of the lens, 
being little sensitive to the lens model used [s ee 
Kochanek's contribution in lMevlan et al. 1 (12001 ], 



2.3. Dynamical Mass 

We call the dynamical mass, mo, the mass esti- 
mated from the observed velocity dispersion. Here 
we are interested in examining the case of a power 



law for the density profile, p — Ar 1 , where A is 
a constant that has to be determined from the 
Jeans equation and the observed velocity disper- 
sion. The mass within the cylinder Ce of Einstein 
radius Re is then 
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mr>(7,/3) 



p(r)dV 



C E 



3 + 7 r(-2) 



? 3+7 



The s pherical Jeans equation (jBinnev fc Tremaine 
19871) can be written as 



1 d(uo?) Pa* ttG 
v dr r dr 3 



A 



.1+7 



'7 



(3) 



where oy is the radial velocity disp ersion , 



19831 



is the luminos ity density profile ( Jaffei 1985 
lHernquistlll990l ). (3 = 1 — of /of is the anisotropy 
parameter of the velocity distribution (<r t is the 
tangential velocity dispersion) , and $ is the grav- 
itational potential produced by the assumed den- 
sity profile. 

Since what is available is the luminosity- 
weighted average velocity dispersion within a 
given aperture, a% , the following constraint is 
necessary for the determination of the constant A: 



(4) 



where the integration volume is an infinite cylinder 
of radius R ap with axis along the l.o.s. 

To simplify, and since there is very little prior 
knowledge on the velocity anisotropy parameter, 
we assume that (3 is a constant. In Appendix lA.il 
we give more details on the solution of the Jean's 
equation ([3]), and in Appendix lA.2l we examine the 
correction due to seeing effects. 

Figure [1] displays the general behavior of the 
dynamical mass as a function of the density pro- 
file slope and velocity anisotropy parameter. We 
examine moipi, = Q) and vtle>(ci = —2, j3). Other 
combinations around these fixed values give quali- 
tatively similar results. The doted line depicts the 
strong lensing mass for this hypothetical system, 
so it is possible to glimpse from the intersection 
of the curves the expected value of the dynamical 
parameters, 7 and (3. The difference of the curves 
calculated using the Jaffe or the Hernquist light 
profiles are just quantitative. The use of a con- 
stant mass to light ratio, v(r) oc p{r) cx r 7 , dis- 
played as the dot-dashed line on the left panel of 
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figure [TJ gives a very distinct and interesting be- 
ha yior for 771.0(7), reproducin g a result obtained 
bv iGuimaraes &, Sodre ( 2007 ). However, a single 
power law is not a realistic approximation for the 
light distribution of the lens-galaxies in the sam- 
ple. 

2.4. Statistical Model 

We want to compare the estimates of lensing 
and dynamical masses taking into account the pos- 
sibility that masses obtained through gravitational 
lensing are affected by some type of contamina- 
tion. Specifically we consider a model where the 
real lensing mass mj, (assumed equal to the dy- 
namical mass) is related to the measured lensing 
mass by 

m L (f e ) = (1 - f c )M L , (5) 

where f c gives the fractional mass contamination. 

Our model, then, has three free parameters: 
f c , 7 and 0. These parameters are determined 
through the Bayes theorem. Since, however, we 
shall assume uniform priors for the parameters 
(see below), the posterior probability of the pa- 
rameters is equivalent to the likelihood. 

To construct a likelihood for the system, which 
has the mission of describing the probabilities of 
the model parameters given the data, we define 
the quantity 



F = 



itid 



m D 
rriL 



(6) 



Note that both and mo refer to the projected 
mass withing the Einstein radius. The likelihood 
for each lens system is then written as 



1 



2na 2 F 



■ exp 



(F t - F 



obs 



(7) 



where F b s = F b s (f c ,j, /3;data) is the measured 
F given the model and the observational data, 
and F t is the expected value for it, which, in 
the desired case where both mj, and mo are esti- 
mates of the same true galaxy mass, corresponds 
to F ( = 0. Note that other quantities may be 
defined to construct the likelihood, for example 
F = m L /m D + m D /m L (F t = 2), F = m L /m D 
(F t = 1), F = m L - m D (F t = 0). The next to 
last definition gives a likelihood that is not sym- 
metrical between and wid, what is not desir- 
able, and the last example has the inconvenience of 



maximizing the likelihood not only in the desired 
region of the parameter space in which mj, ~ rnjj, 
but also in the region where both mass estimates 
are small, what introduces artificial solutions that 
give maximum l.o.s. contamination. However if 
f c = is fixed, then all definitions for F, includ- 
ing the last, give very similar results. 

The variance in F is estimated as being (an 
index i is implied in all quantities) 



(' 



1 d) 



l L u D 



l D u L 



(8) 



D 



where <7{d,l} is the uncertainty in mir),L}- 

The joint likelihood for the whole set of N 
galaxies is then given by 



iV 



£(/ c ,7,/3)=n £ « 



(9) 



which allows us to find the maximum likelihood 
estimator for the free parameters that we denote 
by putting a hat over the parameter (/ c , 7, p), 
and the Bayes estimator (also called the posterior 
mean) defined by (all priors used are flat) 



(P) 



J C(p)pdp 
f£(p)dp ■ 



(10) 



We also calculate the root mean square deviation, 
rms(p) = \J (p 2 ) — (p) 2 , fo r each free p arameter 



and the Bayesian evidence (|Trottall200 



E = J C{p)P{p)dp 



(11) 



where P(p) is the prior probability distribution 
in the parameter space p. We adopt a flat prior 
probability distribution for all the parameters con- 
sidered: f c £ [0,0.5], 7 € [-3,-1] and/3e [-1,1]. 
We sample the parameter space using a grid and 
use it to make the various calculations. 

The main appeal of this approach for model 
comparison is that the Bayesian evidence auto- 
matically implements Occam's razor by penalizing 
more strongly more complex models, those with 
more free parameters. 
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Fig. 1. — Dynamical mass behavior in relation to density profile logarithmic slope, 7, and the velocity 
anisotropy parameter, /3, in units of h Mq. The curves were calculated for an "average" system with 
zi = 0.2, z s = 0.6, (j ap — 250fcms~ 1 , 0e = 1.2", 9 e f = 2.2". Solid lines use Jaffe and dashed use Hcrnquist 
light distribution profiles. The dotted line depicts the strong lensing mass, equation (JXJ) , for the "average" 
system: Ml = 1.34- l0 11 h~ 1 M Q . The dot-dashed line on the left panel was obtained using v(r) oc p(r) oc r 7 . 



3. Results 

We have analyzed the data described in Section 
2.11 with a model with three degrees of freedom, 
one (f c ) corresponding to a possible l.o.s. mass 
contamination that would affect the strong lensing 
estimate of the lens mass and two (7 and 0) cor- 
responding to intrinsic properties of the lens that 
determine its dynamical mass estimate. Two light 
distribution profiles (Jaffe and Hcrnquist) were ex- 
amined. 

Table Q] summarizes the results for the maxi- 
mum likelihood and Bayes estimators for each free 
parameter, as well as for their root mean square 
deviations, which give a measure of their disper- 
sion. For the fractional mass contamination the 
maximum likelihood estimator is zero (no l.o.s. 
contamination), but its mean value is of the order 
of 10%. The density profile logarithmic slope is es- 
timated to be close to an isothermal profile, using 
the Jaffe's light distribution, or somewhat steeper, 
using Hernquist's. In both cases it is found a ~ 5% 
dispersion, w hat is considerably high er than what 
was found bv lKoopmans et all ( 20091 ) (~ 1%), ne- 
glecting l.o.s. contamination and anisotropy. The 



overall effective velocity anisotropy parameter is 
compatib le with the null value (isotropy). We 
note that lKoopmans et al.l ( 20091 ) found a positive 
anisotropy (significantly distinct from isotropy) 
for the same sample, but using for that an inde- 
pendent determination based on scaling relations 
of the density profile logarithmic slope. That is a 
different method from ours, which relies solely on 
the joint strong lensing and dynamical analysis. 

Figure [5] shows the likelihood marginalizations 
over one free parameter. The top row plots use 
Jaffc's profile and the bottom ones, Hernquist's. 
In figure [3] the marginalizations are over two free 
parameters. Qualitatively, the results for both 
profiles are the same. The likelihood distribution 
for the fractional mass contamination indicates a 
peak in f c = with a tail falling to f c ~ 0.5 in the 
3cr confidence region. 

For 7 and j3 the likelihood follows a near Gaus- 
sian behavior, their maximum likelihood and mean 
values are about the same. Jaffe's light distribu- 
tion implies a slightly flatter density profile and 
higher anisotropy compared to Hernquist's. The 
f c marginalized likelihoods (last row of figure [2|) 
show a degeneracy between 7 and (3 that was 



■5 



Table 1: Maximum likelihood parameters, posterior means and root mean squares. 



light profile 


fc 


(fc) 


rms(/ c ) 


7 


(7> 


rms(7) 


a 


(P) 


rms(/3) 


Jaffe 


0. 


0.13 


0.11 


-2.04 


-2.08 


0.12 


0.11 


0.20 


0.21 


Hernquist 


0. 


0.10 


0.09 


-2.18 


-2.13 


0.11 


-0.01 


0.05 


0.17 



already hinted in figure [lj since the curves for 
to_d(7,/3 = 0) and tod (7 = — 2, 0) have mono- 
tonically crescent and decrescent behavior, respec- 
tively. 

If the l.o. s. contamination is neglected, fix- 
ing f c = in our model, there is no significant 
change in the determination of the density profile 
logarithmic slope. However if the velocity distri- 
bution is assumed to be isotropic (/3 = 0), then 
there is a small steepening of the density pro- 
file, 7 = —2.12, and the probability distribution 
dispersion is smaller, rms(7) = 0.06 when using 
Jaffe's light profile, compared with the model with 
three degrees of freedom (Table [TJ . Similarly, for 
Hcrnquist's profile, there is no significant change 
in 7 and rms(7) = 0.07. 

If both fractional contamination and anisotropy 
are fixed, f c = j3 = 0, then 7 = 0.12 for both 
light profiles and rms(7) = 0.019 using Jaffe and 
rms(7) = 0.027 using Hernquist. Therefore, the 
neglect of our ignorance on (3 and, to a lower ex- 
tent, f c , by fixing them equal to zero, implies a 
considerable underestimate of the dispersion in 7 
and possibly the introduction of a systematical er- 
ror. In the central panel of figure [3] it is shown the 
likelihood (narrower curves) for a model with fixed 
f c = (3 = 0, free 7 and Hernquist profile, reproduc 



i ng clo sely what was obtained bv lKoopmans et al 
( 2009h . In the same figure we examine the im- 



pact of the seeing correction (see Appendix IA.2[) 
The two narrow and very close curves differ just in 
that in their calculation one takes into account the 
seeing correction (dot-dashed line) and the other 
(dotted line) does not. The seeing correction is 
negligible. 

To compare models we calculate the ratio of the 
Bayesian evidences, also known as Bayes factor, 



B(M 1 ,M 2 ) 



E(M t 
E{M 2 



(12) 



Table [2] gives a summary of several model com- 
parisons. The Bayes factor between the models 
using Jaffe's or Hernquist's light profile indicates 
that the first is favored over the second. We also 
compare the gain in Bayesian evidence yielded by 
each free parameter separately in relation to a 
model with fixed values (f c — 0, 7 = —2, (3 = 0). 
In general, there is no or weak evidence favoring 
the simpler models (with fixed parameter values) 
over the more complex ones. However, if all pa- 
rameter values are fixed the resulting model with 
no free parameter is strongly (Jaffe) or moderately 
(Hernquist) disfavored over the model with three 
free parameters 

4. Discussion and Conclusions 

Two independent methods, strong gravitational 
lensing and stellar dynamics, were used to deter- 
mine the projected galaxy mass within its Einstein 
radius for a set of 58 galaxies from SLACS. From 
the comparison of the two masses, the lens den- 
sity profile, velocity anisotropy and l.o.s. mass 
contamination were probed. 

Even thought the maximum likelihood esti- 
mation of the l.o.s. mass contamination gives 
a null contribution from secondary matter along 
the l.o.s., the distribution of possible contami- 
nation fraction values is broad, yielding an av- 
erage contamination of the order of 10%. This 
result is compatible with theoretical expecta- 
tions coming from large-scale structure simula- 



tions dPremadi fc Martelll2004lHilbert et al.l2007 



Faure et al.ll2008t IPuchwein fc Hilberd 120091) and 
also with some observational wor ks (Momchev a et al 
20061 : 1 Williams et alJEjoi iTreu et all l2008h that 



Augerl (|2008i rdid not 



which values c an be interpre ted qualitatively using 
Jeffrey's scale (|Trottall2008h . 



do find l.o.s. contamination, 
find an over-density of photometric sources along 
the l.o.s. of a sample of SLACS strong lenses in 
comparison with other SDSS massive early-type 
galaxies. However, it is worth pointing out that 
this is not in contradiction with the existence of 
mass over-densities along the l.o.s. since SLACS 
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7 /J 

Fig. 2. — Marginalized likelihood contours, corresponding to the la, 2a and 3a confidence levels. The 
parameter values of maximum marginalized likelihood are shown as points in the interior of the contours. 
Upper panels use Jaffc's light distribution profile and bottom panels use Hcrnquist's. The marginalization 
is over one of the three free parameters: (3 in the left row, 7 in the center and f c in the right. 




Fig. 3. — Marginalized likelihood over two parameters. The narrower curves at the center panel is for a 
model with fixed f c = /3 = and free 7, using a Hernquist luminosity density profile: the dotted line does 
not take into account the seeing correction and the dot-dashed line does. 
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Tabic 2: Comparison of models. 



model 1: Mi 


model 2: M 2 


C m a,AM 2 ) 


]„ E(Mi) 
E(M 2 ) 


Jeffreys' scale 


J, free (/ c , 7, (3) 


H, free (/ c , 7, f3) 


14 


2.9 


moderate evidence for Mi 


J, free (/ c , 7, /3) 


J, / c = 0, free (7, 0) 


1.0 


-1.1 


weak evidence for M 2 


H, free (/ c , 7, /?) 


H, / c = 0, free (7, 0) 


1.0 


-1.4 


weak evidence for Mi 


J, free (/ c , 7, /3) 


J, 7 = -2, free (/ c , /3) 


1.1 


-1.7 


weak evidence for Mi 


H, free (/ c , 7, /3) 


H, 7 = -2, free (/ c , /?) 


2.4 


-0.69 


inconclusive 


J, free (f c , 7, /3) 


J, (3 = 0, free (/ c , 7) 


1.3 


-0.70 


inconclusive 


H, free (/ c , 7, 0) 


H, /? = 0, free (/ C! 7) 


1.0 


-1.6 


weak evidence for M2 


J, free (/ c , 7, /3) 


J, f c = 0, 7 = -2, /3 = 


4.6-10 7 


11 


strong evidence for Mi 


H, free (/ c , 7, /?) 


H, / c = 0, 7 = —2, /3 = 


1.6-10 4 


3.1 


moderate evidence for Mi 



Note. — J stands for Jaffe's light distribution, H for Hernquist's. Jeffreys' scale according to iTrottal II2008T ). 



lenses were selected to be isolated and excess pho- 
tometric sources would only trace rare high den- 
sity peaks, but not more common diffuse mass 
concentrations. 

The shape of the f c x 7 and f c x (3 degenera- 
cies indicates that the determination of the density 
profile and anisotropy parameter are not very sen- 
sitive to the value of the l.o.s. contamination. This 
can be understood from the functional dependence 
of the dynamical mass with 7 and (3 (figure [J) . A 
small variation in 7 or /? implies a large change in 
m D . 

The joint strong lensing and stellar dynam- 
ics analysis seems to do not strongly constrain 
the anisotropy parameter, its likelihood distribu- 
tion being broad and statistically compatible with 
isotropy on average. Nevertheless, the degeneracy 
between 7 and (3 gives an indetermination that 
correlates a larger anisotropy with a flatter den- 
sity profile, which can also be understood from 
the functional behavior of 77173(7, (3) (figure [l}. An 
increment in 7 can be annulled, in terms of a vari- 
ation in the dynamical mass, by a decrement in /?, 
and vice-versa. 

The inclusion of two degrees of freedom in the 
model with respect to the l.o.s. contamination 
and velocity anisotropy allows us to take into ac- 
count and examine the effect of our prior igno- 
rance on these features of the lens system. The 
most visible effect of f c and (3 on the determi- 
nation of the density profile logarithmic slope 
is the considerable broadening of its likelihood 



distribution (posterior probability), what means 
a more uncertain determination of 7. Within 
this uncertainty, the density profile is statisti- 
cally compatible with an isothermal profile. This 
very particular density profile is also found by 
other authors, apparently as a result of the com- 
pleme ntarity of baryonic a nd dark matter pro- 
files dHamana et all Ejoi iFerreras et al.l f2 005: 



Lintott et al.ll2006HBaltz et al.ll2007t ICzoske et al 
20081) 



We used standard assumptions and approxi- 
mations in our modeling of galaxies and anal- 
ysis. Nevertheless, we can identify several ar- 
eas where further work can be done to refine 
the understanding of these strong lensing systems, 
which can also be seem as caveats to the present 
works in the area. Among them we highlight 
(i) the triaxiality and substructure of lens ha- 
los, whose impor t ance was already suggested by 



los, wnose impor t ance was already suggested by 
iMeneghetti et all (|2005h and lYencho et all (|2006l ) 
in the context of simulations on cluster scales; 
(ii) the correction of the dynamical mass esti- 
mate due to rotational support. It is well known 
that some early-type galaxies ca n have a signifi- 
cant rotational component fe.g.. lEmsellem et al 



20071 ) and two-dimensional kinematics for some 



few SLACS lenses are already becoming available 



(jCzoske et al.l 2008; B arnabe et al.ll2009r ): (iii) the 



brightness distribution, which could be treated 
more realistically with the observed full surface 
luminosities for the individual lenses, instead of 
individualized fits of an universal profile; and (iv) 
the velocity anisotropy parameter, which was as- 
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sumed to be a constant, but that more realistically 
must be a function of radius. However very little 
is known about the velocity anisotropy of early- 
type galaxies, be it observationally, theoretically 
and even from simulations. 

The relaxation of some of our assumptions, with 
the almost inevitable addition of extra free param- 
eters, could prove a fruitful source of investiga- 
tion, however would likely require a larger galaxy 
sample to reduce the likely increased degeneracies 
among the degrees of freedom of the model. As 
we have illustrated, models with a large number 
of parameters may have a higher likelihood but 
lower Bayesian evidence, since the added complex- 
ity must pay its price in a Bayesian sense. 

The authors thank CNPq and FAPESP for fi- 
nancial support and the SLACS and SDSS teams 
for the databases used in this work. 
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Appendix 

A.l. Solving Jean's Equation 

The spherical Jeans equation §3§ can be rewritten, defining x = r/R e f (dimcnsionless) and 

_ 3 + 7 D-(2+7)/-.._2> 



y 



4irGA Re f 



as 



dx x 



(Al) 
(A2) 

x~<*(x + x^-"t* ' (A3) 

where 7* = 1 and = 1/1.8153 (|Hernquistlll"990l ). or 7, = 2 and x* = 1/0.7447 (|Jaffelll983l) . We examine 
both profiles. 

The first-order linear differential equation (|A2|) has as solution 



The luminosity distribution is well approximated by the profiles ([Koopmans et al.ll200(: 

v(x) cx 



y(x) = x-^ 



C 



iy(x)x 1+ ^ +2f3 dx 



(A4) 



where C is an arbitrary constant and the most evident boundary condition is y(x — > 00) = 0. The analytical 
solution for the integral in (|A4[) with (|A3|) for v(r) has the hypergeometric function 2 -Pi, which has a 
computationally demanding solution. Therefore we solve (|A2|) using a fourth order Runge Kutha algorithm, 
starting at y{x — 1000) = and evolving y{x) down to x = 0.01. 
The dynamical mass ([2]) within the Einstein radius Re is then 



Re\ 2+1 U, vdV 



R ef J Icu y dv 



(A5) 



and its error is estimated from the observational errors on a ap and Re through error propagation. 
A. 2. Seeing 

We model the effect of the seeing through a Gaussian smoothing of the galaxy projected luminosity. 
Therefore, the observed surface brightness profile is related to an intrinsic (no seeing) profile by 



Iobs{8) 



i(e')i 



(A6) 



where Io is the modified Bessel function of first kind, and a 2 is the Gaussian seeing variance. We use 
a s = 0.64", which corresponds to a FWHM of 1.5". 

The seeing correction of the average velocity dispersion within the observational aperture is then given by 



Io ap o$(R)I(R)RdR f* ap I obs (R)RdR 
[o*(R)I(Rj\ obs RdR J*"* I(R)RdR 



a P> obs 



(A7) 



where a 2 i s the projected velocity dispersion profile, [ a 2 (R)I(R)] ^ is defined in an analogous way to I bs (R) 
in Eq. (|A6P and the projection is calculated through ( Binnev fc Tremaine 1119871 ) 



Vr 2 - R 2 



zdr. 



(A8) 
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